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TIME VARIATION IN G24.78+0.08 Al: EVIDENCE FOR AN 
ACCRETING HYPERCOMPACT H ii REGION? 

Roberto Galvan-Madrid^'^, Luis F. Rodriguez^, Paul T. P. Ho^'^, and Eric Keto^ 

ABSTRACT 

Over a timescale of a few years, an observed change in the optically thick 
radio continuum flux can indicate whether an unresolved H II region around a 
newly formed massive star is changing in size. In this Letter we report on a study 
of archival VLA observations of the hypercompact H II region G24.78+0.08 Al 
that shows a decrease of ~ 45 % in the 6-cm flux over a five year period. Such 
a decrease indicates a contraction of ~ 25 % in the ionized radius and could be 
caused by an increase in the ionized gas density if the size of the H II region 
is determined by a balance between photoionization and recombination. This 
finding is not compatible with continuous expansion of the H II region after the 
end of accretion onto the ionizing star, but is consistent with the hypothesis of 
gravitational trapping and ionized accretion flows if the mass-accretion rate is 
not steady. 

Subject headings: H II regions — ISM: individual (G24.78+0.08) — stars: for- 
mation 

1. Introduction 

The formation of massive stars (M > 8 Mq) by accretion presents a number of theoreti- 
cal difficulties; among them, that once a star attains a sufficient mass its surface temperature 
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is high enough to produce a small H II region (see the reviews on ultracompact and hyper- 
compact H II regions by Churchwell 2002; Kurtz 2005; and Hoare et al. 2007). The thermal 
pressure differential between the hot (~ 10^ K) ionized gas and the cold (~ 100 K) molecu- 
lar gas can potentially reverse the accretion flow of molecular gas and prevent the star from 
ever reachi ng a higher mass. However, recent models that include the effects of gravity (e.g., 
Ketdl2007l ) have shown that if the H II region is small enough (ionized diameter < 1000 AU 
for total stellar mass > 350 M©), the gravitational attraction of the star(s) dominates the 
thermal pressure and the molecular accretion flow can cross the ionization boundary and 
proceed toward the star as an ionized accretion flow within the H II region. In this stage the 
H II region is said to be gravitationally trapped by the star. 

Previous studies have demonstrated a number of observational techniques bearing on the 
evolution of small H II regions around newly formed stars. Accretion flows onto and through 
H II regions can be directly observed by ma pping molecular and radio recombinati on lines 
(RRLs) at very high angular resolution (e.g.. ISollins et al.ll2005l : iKeto fc Woodll2006l ). If the 
H II region is too small to be spatially resolved, the frequency dependence of the velocities 
and widths of RRLs ca n be used to infe r steep density gradients and supersonic velocities 
within the H II regions (IKeto et al.ll2007l ). as expected from the presence of ionized accretion 
flows or bipolar outflows. Radio continuum observations at very high angular resolution, 
made at two epochs some years apart, can also be used to directly observe changes in size of 
small H II regions (Franco- Hernandez & Rodriguez 2004). Changes in size indicate whether 
the evolution of the H II region is consistent with pressure-driven expansion or gravitational 
trapping. 

In this Letter we demonstrate yet another technique - how a comparison of radio con- 
tinuum observations at different epochs can be used to infer size changes in H II regions 
even if the observations do not spatially resolve them. At optically thick frequencies the 
radio continuum flux depends to first approximation only on the size of an H II region and 
is independent of its internal density structure. Therefore, even if the size is not known, 
a change in the flux over time still indicates a change in the size of the H II region. This 
technique is particularly valuable because it can be used on the smallest and youngest H II 
regions, and because lower angular resolution radio continuum observations generally require 
less observing time than spectral line and high angular resolution observations. 

For the study presented in this Letter we selected the hypercompact (H C) H II region 



G24. 78+0.08 Al which lies at the center of a massive molecular accretion flow ( Beltran et al. 



2004 . 12OO6I ) and also has multi-epoch 6 cm radio continuum observ ations in the 



G24. 78-1-0. 08 wa s detected in the cen timetric (cm) continuum by iBecker et al.l (119941 ). and 



J A archive. 



later resolved by ICodella et al.l (119971 ) into a compact (A) and an extended (B) component. 
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A millimeter (mm) interfero metric study (IBeltran et al.ll2004l ) revealed the presence of two 
massive rotating toroids centered in respective dust cores (Al and A2). The compact cm 
emission co mes from the r nm c omponent Al (hereafter G24 Al), and has recently been 
resolved by iBeltran et al.l ( 20071). If G24 Al is ionized by a single star, its spectral type 
should b e earlier than 09 (ICodella et al.l 119971 ) . Also, G24 Al likely powers a massive CO 
outflow JPuruva et al.ll2002h . 

The infalling and rotating molecular gas and bipolar outflow all suggest ongoing accre- 



tion. Howe ver, based on proper mo tions of H2O masers around the H II region, IBeltran et al. 



( I2OO7I ) and iMoscadelli et al.l (120071 ) proposed that at the present time, the H II region is ex- 
panding into the accretion flow. The suggested timescale for the expansion is short enough 
that we should be able to detect a corresponding increase in the optically thick radio con- 
tinuum flux within a few years. 



2. Observations 

We searched the VLA0 archive for multi-epoch observations centered in the G24 Al 
region at optically thick frequencies (below 23 GHz for G24 Al). Since in the optically thick 
part of the spectrum the flux density of the source scales as the angular size squared, flux 
density variations corresponding to size variations can be detected even in observations of 
modest angular resolution. We chose three data sets of 6-cm observations in the C configu- 
ration, two from 1984 and one from 1989 (see Table 1). 

The observations were made in both circular polarizations with an effective bandwidth 
of 100 MHz. The amplitude scale was derived from observations of the absolute amplitude 
calibrator 3C286. This scale was transferred to the phase calibrator and then to the source. 
We estimate an error not greater than 10 % for the flux densities of the sources. 

We edited and calibrated each epoch separately following the standard VLA procedures 
using the reduction software AIPS. Precession to J2000 coordinates was performed running 
the task UVPIX in the {u, v) data. After self- calibration, we made CLEANed images with 
uniform weighting and cutting the short spacings (up to 10 KA) to minimize the presence of 
extended emission at scales larger than ~ 20". 

Before subtraction, we made the images as similar as possible. We restored the CLEAN 
components with an identical Gaussian beam HPBW 4''79 x 3''38, PA= —15°. We applied 
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primary beam corrections and aligned the maps. No significant differences were found in 
the difference image between the 1984 May 11 and 14 images, as expected for such a small 
time baseline. We therefore averaged the two 1984 epochs into a single data set. In sub- 
tracting the final maps we allowed the 1984 image to have small (~ 0.1") shifts in position 
as well as a scaling of ~ 10 % in amplitude. This was done in order to minimize the rms 
residuals of the difference image in th e region of interest. A similar procedure was used by 
Franco-Hernandez fc Rodriguezl (120041 ) to detect a variation in the lobes of the bipolar UC H 
II region NGC 7538 IRSl. The individual maps, as well as the final difference image between 
1984 and 1989 are shown in Figure 1. 



Discussion 



3.1. The Expected Variation Trend 
3.1.1. If G24 Al Is Expanding 



Beltran et al.l (120071 ) suggested that the 7 mm and 1.3 cm continuum morphologies are 
consistent with limb-brighten ing from a thin, ion ized-shell structure. Based on H2O maser 
proper motions (see Fig. 3 of iBeltran et al.l 120071 . or Fig. 4 of Moscadelli et al. 2007) they 
also suggested an expansion speed of ~ 40 km s~^. The increase in flux corresponding to 
the increase in size due to expansion ought to be detectable in a few years. 



3.1.2. If G24 Al Is Accreting 

If G24 Al is the ionized inner portion of the star forming accretion flow, the long- 
term growth of the H II region due to the increasing ionizing flux of the star should be 
imperceptible. However, the H II region could change in size over an observable timescale if 
the gas density in the accretion flow is time variable. Because the mass of ionized gas within 
the H II region is very small compared to the mass of the accretion flow, even a small change 
in the flow density could affect the size of the H II region. 

For example, we can obtain a lower limit for th e mass of G24 Al ass u ming that it is 
spherical and homogeneous. From the equations of iMezger fc Henderson! (119671 ) we have 
that: 
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where the flux density has been measured at a frequency v in which the H II region is 
optically thin, is the electron temperature, D is the distance to the region, and 6. is it s 
FWHP. Taking = 101 mJy, D = 7.7 kpc, and 9s = O'.'U at 7 mm JSeltran et al.lboOTh . 
and assuming = 10^ K, the ionized mass in G24 Al would be ~ 5 x 10~^Mq. A higher 
limit to the mass can be obtained considering that HC H II regions shou ld have den sity 
gradients (n oc r°, with a = —1.5 to —2.5) rather than being homogeneous faeto 2007 ). In 



this case we obtain ionized masses between 1 x 10 ^ Mq and 3 x 10 ^ M^, a factor of 2 — 6 



'05 



higher than under the assumption of homogeneity, but still very small when compared to 
that of the surrounding molecular material or the ionizing star itself. 



3.2. The Observed Variation 

Figure 1 shows the individual maps of the 1984 May 11 + 14 (1984.36) and 1989 June 
23 (1 989.48) epochs, as w ell as the difference image. Sources Al and B (the former labeled 



A by lCodella et al.l 119971 ) are unresolved {6s < 2"), and the total emission is dominated by 
component B. However, the flux decrease is centered at the position of Al. We performed 
Gaussian fits to both components using the task JMFIT in AIPS, and the results are sum- 
marized in Table 2. We have checked the reliability of the fits to the image by making direct 
fits to the (u,v) data. The values obtained from both techniques are entirely consistent but 
suggest that the errors given by the tasks used in the image (JMFIT) and (u,v) (UVFIT) 
fittings are underestimated by a factor of 2. The errors given in Table 2 have been corrected 
by this factor. 

The decrease in the flux density of G24 Al between 1989.48 and 1984.36 is 5.1 ± 1.1 
mJy, or 45 ± 10 %. Component B shows no evidence of time variability, as expected for a 
more evolved H II re gion not associated with signs of current star-forming activity such as 



H2O or OH masers (jCodella et al.l 119971 ). We set an upper limit of 2 % to the circularly- 
polarized emission of the sources, which indicates that we are not dealing with variable 
gyrosynchrotron emission from an active stellar magnetosphere. 



3.3. Is G24 Al Accreting? 



The ~ 45 % flux decrease (i.e., a contraction of ~ 25 % in the ionized radius) we have 
detected at 6 cm toward G24 Al is not consistent with the hypothesis that this H II region 
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is expanding rapidly into the molecul ar accretion flow. A ssuming a radius R ~ 500 AU and 
an expansion velocity v ~ 40 km s"^ flBeltran et al.l 120071 ). a flux increase of ~ 20 % should 
have been observed between the compared epochs. 

We attribute the contraction of the H II region to an increase in its density produced by 
the enhancement of accretion, either caused by an isotropic increment in the mass-accretion 
rate or by the sudden accretion of a localized clump in the neutral inflow (for e xample, clumps 
have been observed in the accretion flow onto the UC H II region GlO.6-0.4: ISoUins and Ho 
20051 : iKeto fc Woodll2006l ). Our data do not allow us to distinguish among these two possi- 
bilities. 

The additional mass required to increase the density can be estimated. The ionized 
radius scales roughly with the density n as Ts oc ?t,~^/^ (if the radius of the H II region is set 
by a balance between photoionization and recombination) and at optically thick frequencies 
the flux density scales as Si, cx r^. The ionized mass within G24 Al is ~ 1 x 10~^ Mq 
( §3.1.2p . Therefore, the sudden accretion of ~ 5 x 10~^ Mq into the H II region would suffice 
to explain the observed 45 % flu x decrease. The mo lecular accretion flow around this H II 
region has a mass of ~ 130 Mq (IBeltran et al.ll2004j ). Thus, the required variation, 50 % of 
the mass of the H II region, is only 0.001 % of the total mass of the accretion flow. 

Finally, we shall mention that even when the observations here reported are not con- 
sistent with a simple, continuous expansion of G24 Al, they do not rule out other more 
complex, non-steady evolutionary scenarios. 



Conclusions 



Our analysis of archival VLA observations of the H II region G24.78+0.08 Al indicates 
a contraction of its radius between 1984.36 and 1989.48. This flnding is consistent with the 
hypothesis that this HC H II is the inn er, ionized pa r t of the lar ger scale accretion flow 
seen in the molecular line observations of Beltran et al. (j2004 . 2006 ). Future high- resolution 
RRL observations could unambiguously resolve the ionized-gas kinematics and conflrm this 
hypothesis. 
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Table 1. Observational Parameters 



Epoch 


Phase Center'' 


Amplitude 


Phase 


Bootstrapped Flux 


Beam 




a(J2000) 


(5(J2000) 


Calibrator 


Calibrator 


Density (Jy) 


(arcsec x arcsec; deg) 


1984 May 11 


18 36 12.145 


-07 11 28.17 


30286 


1743 - 038 


2.414 ± 0.005 


5.38 X 3.44; -29 


1984 May 14 


18 36 12.145 


-07 11 28.17 


3C286 


1743 - 038 


2.464 ± 0.007 


4.78 X 3.34; +03 


1989 Jun 23 


18 36 10.682 


-07 11 19.87 


3C286 


1834 - 126 


0.153 ± 0.001 


4.50 X 3.31; -10 



*Units of right ascension axe hours, minutes, and seconds. Units of declination are degrees, arcminutes, and arcseconds. 
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Table 2. Gaussian Fit Results 



Epoch 


Component 


Position^"'' 


Flux Density 






a(J2000) 


6(J2000) 


(mjy) 


1984.36 


Al 


18 36 12.545 


-07 12 10.87 


11.4 ±0.8 


1984.36 


B 


18 36 12.668 


-07 12 15.37 


31.4 ± 1.0 


1989.48 


Al 


18 36 12.. 54.5 


-07 12 10.87 


06.3 ± 0.8 


l<),s<).is 


B 


IS 12.()()S 


-07 12 l."').:-!7 


:-l2..s = 0,8 



■^Units of right ascension are hours, minutes, and seconds. Units of dech- 
nation are degrees, arcminutes, and arcseconds. 

''For the 1989.48 epoch, the centers of the Gaussians were fixed to those 
obtained for the 1984.36 epoch. 
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Fig. 1.— VLA images of G24.78+0.08 for 1984.36 (top), 1989.48 (middle), and the difference 
of 1989.48 - 1984.36 [bottom). The contours are -10, -8, -6, -5, -4, 4, 5, 6, 8, 10, 12, 15, 20, 
30, 40, and 60 times 0.57 mJy beam~^. The half power contour of the synthesized beam 
(4''79 X 3''38 with a position angle of —15°) is shown in the bottom left corner of the images. 
The crosses indicate the positions of the components Al and B from our Gaussian fits to the 
1984.36 image. The negative residuals observed in the difference image indicate a decrease 
of ~ 45 % in the flux density of component Al. 



